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ABSTRACT 

We present new, intermediate resolution spectra (~ 4 A) of the compact radio source 
PKS 1345+12 (4C 12.50, z = 0.122) with large spectral coverage (- 4500 A). Our 
spectra clearly show extended line emission up to ^ 20 kpc from the nucleus. This is 
consistent with the asymmetric halo of diffuse emission observed in optical and infra- 
red images. At the position of the nucleus we observe complex emission line profiles. 
Gaussian fits to the [O III] emission lines require 3 components (narrow, intermediate 
and broad), the broadest of which has FWHM ^ 2000 km s^^ and is blue shifted 
by up to ~ 2000 km s~^ with respect to the halo of the galaxy and HI absorption. 
We interpret this as material in outflow. We find evidence for high reddening and 
measure E(B-V) > 0.92 for the broadest, most kinematically disturbed component. 
This corresponds to an actual H/3 flux 130 times brighter than that observed. From 
our model for [S II]AA6716,6731 we estimate electron densities of Ug <150 cm~'^, Ug 
>5300 cm~^ and Ug >4200 cm~^ for the regions emitting the narrow, intermediate 
and broad components respectively. We calculate a total mass of line emitting gas 
of Mgas < 10^ Mq. Not all emission line profiles can be reproduced by the same 
model with [O I]AA6300,6363 and [S II]AA6716,6731 requiring separate, unique mod- 
els. We argue that PKS 1345-1-12 is a young radio source whose nuclear regions are 
enshrouded in a dense cocoon of gas and dust. The radio jets are expanding through 
this cocoon, sweeping material out of the nuclear regions. Emission originates from 
three kinematically distinct regions though gradients (e.g. in density, ionisation po- 
tential, acceleration etc) must exist across the regions responsible for the emission of 
the intermediate and broad components. 

Key words: ISM: jets and outflows - ISM: kinematics and dynamics - galaxies: 
active - galaxies: ISM - galaxies: kinematics and dynamics - galaxies: individual: PKS 
1345+12 (4C12.50) 



1 INTRODUCTION 



Gigahertz Peaked Spectrum radio sources (GPS: D < 1 
kpc) and the more extended Compact Steep Spectrum radio 
sources (CSS: D < 15 kpc) make up a significant fraction 
(~ 40%) of the radio source population, but their nature 
is not fully understood (see e.g. O'Dea 1998 and references 
therein). Are these compact radio sources linked to the large 
scale radio galaxies by some evolutionary sequence, or are 
they indeed a separate class of object? 

There are two main theories to describe compact ra- 
dio sources: the frustration scenario (van Breugel 1984) and 
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the youth scenario (Fanti et al. 1995). In the frustration sce- 
nario, the radio source is enshrouded by a cocoon of gas and 
dust so dense that the radio jets cannot escape the nuclear 
regions and the radio source is confined and frustrated for its 
entire lifetime. Alternatively, the youth scenario attributes 
compactness of the radio source to evolutionary stage - if 
we observe a radio source when it is relatively young, the 
radio jets will not have had the time to expand much and 
will still be relatively small. In this scenario, the GPS and 
CSS sources represent the earliest stages in the evolution 
of radio galaxies, eventually expanding into the large scale 
doubles. This is currently the preferred theory and is sup- 
ported by estimates of dynamical ages of GPS sources of 
tdyu ~ 10^ - 10^ years (e.g. Owsianik, Conway & Polatidis 
1998; Tschager et al. 2000) and radio spectral ages of the 
larger CSS sources of t^p < 10* years (Murgia et al. 1999). It 
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Table 1. Log of observations. The red and blue arms were exposed simultaneously using a dichroic at 6100 A. 
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has also been suggested that the host galaxies pass through 
multiple, short-lived (trs < 10* years) activity phases re- 
sponsible for the extended, diffuse emission commonly found 
around GPS sources (Rcadhcad ct al. 1994; Readhead et al. 
1996), although there arc no strong arguments either for or 
against this theory (O'Dca 1998). 

There is much evidence to suggest that, in a significant 
fraction of the radio source population, the radio sources are 
triggered by merger events involving two or more galaxies, 
of which at least one is gas rich (Heckman et al. 1986). Many 
compact radio sources exhibit features attributed to mergers 
such cis double nuclei, tidal tails, arcs of emission and dis- 
torted isophotes, suggesting that these sources are possibly 
observed relatively shortly after a merger event (Stanghellini 
et al. 1993), before the system has recovered and settled 
down. However, alternative activity triggering mechanisms 
have been suggested such as intracluster cooling flows (Bre- 
mer, Fabian & Crawford 1997). Whilst mergers may be the 
most likely trigger of the activity, it is clear that all activ- 
ity inducing processes involve the injection of substantial 
amounts of gas and dust into the nuclear regions - initially, 
at least, the radio source will exist in a dense and dusty 
environment. Subsequent evolution will be strongly depen- 
dent on both the quantity and distribution of the circum- 
nuclear material. If the ISM forms an equatorial disk, the 
radio jets will undergo little interaction with the ISM. On 
the other hand, if the ISM forms an enshrouding cocoon, 
strong interactions are expected between the ISM and the 
radio jets (Bicknell & Dopita 1997), quasar-induced winds 
(Balsara & Krolick 1993) or any starburst-driven superwinds 
(Heckman, Armus & Miley 1990) in the early stages of ra- 
dio source evolution. Indeed, depending on the amount of 
material in the cocoon, it may be necessary to hollow out 
bi-polar cavities in the cocoon in order to allow the radio 
source to expand freely through the halo of the host gala^xy. 
Therefore, it seems likely that outflows of the ISM will be 
important in the early stages of radio source evolution. 

Much of the work on compact radio sources has concen- 
trated on the radio wavelength region, but relatively little 
attention has been paid to other wavebands. High quality 
optical spectra have the potential to provide important clues 
on the nature of compact radio sources and their early evo- 
lution. However, to date, most spectroscopic studies have re- 
lied on poor quality, low resolution data (e.g. Gelderman & 
Whittle 1994). Recently, the promise of optical spectroscopic 
studies was demonstrated by observations of the compact, 
flat spectrum radio source PKS 1549-79 (Tadhunter et al. 
2001). Low dispersion spectra of this radio source revealed 




Figure 1. Slit positions for the new spectroscopic observations 
of PKS 1345H-12 superimposed on an O III image taken with 
HST WFPC2 from Axon et al. (2000). Radio emission is confined 

entirely within the western nucleus (at the centre of the image) 
on a scale of < 0.1 arcsec (~ 240 pc; Stanghellini et al. 1993, 
Stanghellini et al. 1997). The second nucleus lies ~ 1.8 arcsec 
(4.3 kpc) to the east and is included in our spectra taken at PA 
104. 

broad forbidden line proflles (FWHM ~ 1350 km s"^) blue- 
shifted by ^ 600 km s~^ with respect to the galaxy rest 
frame, believed to trace material in outflow. Combined with 
its compactness, Tadhunter et al. (2001) suggest PKS 1549- 
79 is a young radio galaxy whose radio jets, aligned close to 
the line of sight, are sweeping material out of the nuclear 
regions. 



In this paper we present new, high quality observations 
of the GPS source PKS 1345-^12, taken with ISIS, the dual- 
arm spectrograph on the 4.2m William Herschel Telescope 
(WHT) on La Palma. The high quality of the data for this 
source allows us to investigate the putative circum-nuclcar 
outflows in greater depth than was possible for PKS 1549- 
79. The preliminary results of this study were reported in 
Holt, Tadhunter & Morganti (2003). ISIS is ideal for studies 
of this type as it allows spectra to be taken with intermedi- 
ate resolution without compromising the spectral range. By 
using a dichroic, the red and blue regions of the spectrum 
were imaged simultaneously with CCDs optimised for each 
range. We investigate the kinematics of the emission line gas 
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in the nucleus with respect to the galaxy halo to search for 
outflows - a possible impact of the nuclear activity on the 
ambient ISM. We also investigate the physical properties in 
the nuclear regions (e.g. density, gas mass) and outline a 
model for compact radio sources using PKS 1345+12 as an 
example of a radio source in the early stages of radio source 
evolution. 



2 PREVIOUS OBSERVATIONS OF PKS 
1345+12 



The bright, compact radio source PKS 1345+12 is opti- 
cally identified with a 17th magnitude elliptical galaxy (e.g. 
Clarke, Bolton & Shimmins 1966; Burbidgc 1970; Gilmorc 
& Shaw 1986). PKS 1345+12 has two nuclei separated by 
~ 1.8 arcsec or 4.3 kpc^. The western nucleus has a Seyfert 
2/NLRG type spectrum (Gilmorc & Shaw 1986). The op- 
tical identification of the radio source is controversial but 
the latest astrometry by Stanghellini et al. (1993) associates 
the radio source with the western, highly reddened nucleus 
(Evans et al. 1999, Smith & Heckman 1989, Stanghellini 
et al. 1993). 

Both nuclei are surrounded by an extended, asymmet- 
ric halo detected in optical (Heckman et al. 1986; Gilmore 
& Shaw 1986) and near infra-red (Evans et al. 1999) images. 
There is clear evidence for extended, distorted morphology 
including a strongly curved tidal tail terminating ~ 25 arc- 
sec (^ 59 kpc) NW of the nucleus (Heckman et al. 1986, 
Smith & Heckman 1989). Gilmore & Shaw (1986) estimate 
a projected linear size of 110 kpc (where h is the Hubble 
constant with units of 100 km s^ Mpc~^). 

There is much evidence to suggest that the nuclear re- 
gions of PKS 1345+12 contain large amounts of gas and 
dust. Deep, multiple component HI 21cm absorption over a 
velocity range > 700 km s~^ with column density NH//Tspf„ 
= 6.1 X 10^* atoms cm~^ has been detected by Mirabel 
(1989) and more recently over a velocity range > 1000 km 
s~^ by Morganti et al. (2002). The same velocity components 
are also detected in CO (1— >0) emission (Mirabel, Sanders 
& Kazes 1989, Evans et al. 1999) and the shared CO halo 
indicates the galaxy is rich in highly concentrated dust and 
molecular gas (3.3 x 10^° M©; pgas > 2000 Mq pc"^; Evans 
et al. 1999). 

Both the distorted optical morphology and the presence 

of a rich ISM suggest that PKS 1345+12 has been involved 
in a merger of two galaxies in its recent history (Gilmore 
& Shaw 1986, Heckman et al. 1986). At least one of the 
galaxies must have been gas rich (Stanghellini ct al. 1993, 
Gilmore & Shaw 1986). Indeed, the close proximity of the 
two nuclei may indicate that the merger event began so re- 
cently that it may still be ongoing (Gilmore & Shaw 1986, 
Heckman et al. 1986, Hutchings 1987, Baum, Heckman & 
van Breugel 1988, Evans et al. 1999). PKS 1345+12 also 
has a significant young stellar population (Robinson et al. 
2003, in preparation) which may bo the result of a merger 
induced starburst. Further evidence for the prodigious star 

^ Ho = 75 km s~^, go = 0.0 assumed throughout, resulting in a 
scale of 2.37 kpc arcsec"^ at 2 = 0.122. 



formation in PKS 1345+12 is provided by the detection of a 
substantial far-IR excess (Lir — 1.7 x 10^^ Lq, Evans et al. 

1999) which leads to the classification of this source as an 
ultra luminous infra-red galaxy (ULIRG). 

The main radio source associated with PKS 1345+12 
is confined to a region < 0.15 arcsec (-^ 350 pc) in size 
coincident with the western nucleus (Stanghellini et al. 
1993, Stanghellini et al. 1997). PKS 1345+12 has a typical 
GPS radio spectrum, peaking at ~ 0.6 GHz with spectral 
indcx^ -0.7 and +0.9 above and below the peak respectively 
(Stanghellini et al. 1998) and a sharp spectral cut-ofi^ near 
400 MHz (Lister, Kellermann & PauUny-Toth 2002a). 

PKS 1345+12 has distorted, 'triple' radio morphology 
(i.e. a core and two jets are visible) aligned along PA 160 
with a jet extending 0.04 arcsec (95 pc) to the SE before 
bending and expanding into a diffuse lobe. Weak radio emis- 
sion is detected up to 0.03 arcsec (70 pc) NW of the bright 
knot towards the northen limit of the jet, tentatively identi- 
fied as the core (Stanghellini et al. 1997). Stanghellini et al. 
(2001) also report very weak, extended emission on axc- 
second scales at 1.4 GHz resembling the parsec scale mor- 
phology. Such extended weak emission may indicate PKS 
1345+12 has experienced other activity phases in the past 
(e.g. Stanghellini ct al. 2002, private communication). 

The radio morphology of PKS 1345+12 resembles only a 
few GPS sources (see samples listed in O'Dea 1998), though 
this may be a redshift issue - PKS 1345+12 is one of the clos- 
est GPS sources (z — 0.122) whilst most are at much higher 
redshift (z > 0.5). The S-shaped jet morphology is only de- 
tected in 5-10 % of GPS sources and is attributed to recent 
AGN outflows where the black hole spin is still precessing 
(Stanghellini et al. 2001, Lister et al. 2002b). Lister, Keller- 
mann & Pauliny-Toth (2002a) detect superluminal motion 
(1.2c) in the southern jet. 

Consistent with its status as a powerful radio source 
{Pr,GHz = lO^*' W Hz"^; O'Dea 1998 and references therein), 
the nucleus of PKS 1345+12 emits strong emission lines. 
However, although classified as an NLRG by Grandi (1977), 
it has very broad forbidden emission lines with complex 
structures and strong blue asymmetries. 

(i) PKS 1345+12 has two distinct redshift systems: one 
for the high iornsation lines (e.g. [O III], [Nc III]), another 
redshifted by +0.0008 producing the hydrogen Balmer and 
low ionisation lines (e.g. [O I], [O II], [S II] and [N II]). 

(ii) The emission lines are much broader than in normal 
NLRGs and exhibit distinct blue wings with the low ioni- 
sation lines (FWHM 1200 km s^^) having different line 
profiles to the high ionisation forbidden lines and hydrogen 
Balmer lines (FWHM ~ 1600 km s"^). 

In addition to the high ionisation optical forbidden line 

emission, further evidence that PKS 1345+12 contains a 
powerful quasar nucleus is provided by the detection of a 
broad (FWHM ~ 2600 km s"^) component to Paa (VeiUeux, 
Sanders & Kim 1997), the detection of a point source compo- 
nent in high resolution near-IR images (Evans et al. 1999), 
and the detection of a nuclear X-ray source (O'Dea et al. 

2000) . 
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Figure 2. Composite spectrum of tlie nucleus (from PA 160; see sections 3 and 4) at rest of PKS 1345+12. 



3 NEW SPECTROSCOPIC OBSERVATIONS 
OF PKS 1345+12 

We present new observations taken in May 2001 with ISIS, 
the dual-arm spectrograph using the 4.2m William Hcrscliel 
Telescope (WHT) on La Palma. In the red, the R316R grat- 
ing was used with the TEK4 CCD, and in the blue, the 
R300B grating was used with the EEV12 CCD. A dichroic 
at 6100A was used to obtain spectra with wavelength range 
6212-7720A in the red and 3275-6813A in the blue simulta- 
neously. A summary of observations is presented in Table 1. 
To reduce the effects of differential refraction, all exposures 
were taken when PKS 1345+12 was at low airmass (sec z < 
0.1) or along the parallactic angle. 

Spectra were taken along three slit positions - PA 104, 
PA 160 and PA 230 (see Figure 1) with a 1.3 arcsec slit. PA 
104 includes the second rmclcus 1.8 arcesc (4.3 kpc) to the 
east of the active nucleus. PA 160 lies along the radio axis 
and cuts through the arc of emission 1.2 arcsec (2.8 kpc) to 
the north of the nucleus. The final position, PA 230, was the 
parallactic angle of the source at the time of observation. All 
angular scales are taken from Axon et al. (2000). 

The data were reduced in the usual way (bias subtrac- 
tion, flat fielding, cosmic ray removal, wavclcugth calibra- 
tion, fiux calibration) using the standard packages in IRAF. 
The two-dimensional spectra were also corrected for spatial 
distortions of the CCD. To reduce wavelength calibration er- 



rors due to flexure of the telescope and instrument, separate 
arcs were taken at each position on the sky. The final wave- 
length calibration accuracy, calculated using the standard 
error on the mean deviation of the night sky emission lines 
fi'om published values (Osterbrock, Pulbright & Keane 1996) 
was O.O59A and O.II2A in the red and blue respectively (PA 
160). The spectral resolution, calculated using the widths of 
the night sky emission lines was 3.66 ± 0.09A in the red and 
4.54 ± 0.10 A in the blue. 

Comparison of several spectrophotometric standard 
stars taken with a wide slit (5 arcsec) throughout the run 
gave a relative flux calibration accurate to ± 5 per cent. This 
accuracy was conflrmed by good matcliiug in flux between 
the red and blue spectra. Further observations of standard 
stars with a narrow slit, matched to the slit width used to 
observe PKS 1345+12, were used to correct for atmospheric 
absorption features (e.g. A and B bands at ~ 7600A and ~ 
6800 A respectively). 

The main aperture used was the nuclear aperture - 2.16 
arcsec wide, centred on the nuclear continuum emission 

The spectra were extracted and analysed using the Star- 
link packages FIGARO and DIPSO. 
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Figure 3. Spatial variations of the [O II] radial velocities of the extended gaseous halo along all three PAs: PA 160 
(open diamonds), PA 104 (open squares) and PA 230 (asterisks). The open diamonds, squares and asterisks represent 
the spatially extended comijoiicnt of [O II] (the error bars arc small and are incorporated in the points). Overplotted 
is the radial velocity of the deep HI 21cm absorption (filled circle; Mirabel 1989) and the three components of [O 
III] in the nucleus (open triangles; see section 4.2). Note, the orientations of the three PAs are as follows: i) PA 
160: positive kpc = south, negative kpc = north; ii) PA 104: positive kpc = east, negative kpc = west; iii) PA 230: 
positive kpc = south-west, negative kpc = north-east. The vertical dotted line represents the centroid of the nuclear 
aperture. The horizontal dotted line represents the velocity at the position of the centroid. 



4 RESULTS 

In the following sections we will describe the results on the 
kinematics and physical conditions derived from the data. 
Discussion of the emission line ratios and ionisation state of 
the gas will be addressed in a future paper. 

4.1 Kinematics of the extended gaseous halo 

A key aspect of this project is to investigate the impact of 
the activity on the circumnuclear material. It is therefore 
necessary to establish the exact redshift of the galaxy rest 
frame. The heliocentric redshift, calculated using the narrow 



component of [O III] in the nuclear aperture (see section 4.2) 
is z = 0.12351 ± 0.00008. 

The [O II]AA3726,3728 emission lines are highly ex- 
tended, detected up to ~ 12 kpc to the N, E, W and SW; 
up to ~ 15 kpc to the NE and up to ~ 20 kpc to the S. Ra- 
dial velocity profiles were constructed by fitting up to three 
Gaussians (narrow, intermediate and broad components; see 
section 4.2) to the doublet along the entire spatial range and 
are presented in Figure 3. For simplicity, [O II] was treated 
as one line rather than a doublet, hence the weighted mean 
of the rest wavelengths was used as the rest wavelength of 
the doublet when calculating the velocity shifts. 

There are no clear patterns in the spatial variations of 
the [O II] radial velocity - there is no steep gradient across 
the nucleus and the overall shape does not resemble the 
signature 's' shape expected for circular motions in a disk 
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Table 2. Line data for all components of emission lines (where the emission line is a doublet, data for the brightest 
line is tabulated) measured in the nucleus of PKS 1345+12 (see section 3). Columns are: (a) emission line; (6) 
emission line components where n, i and b are the narrow, intermediate and broad components respectively; 
(c) the quoted rest wavelength of the line in angstroms; (d) velocity shift of each component with respect to 

the narrow component for a particular emission line (km s^^); (e) uncertainty on the velocity shift (km s^^); 
(/) velocity width (FWHM, corrected for instrumental width) of the emission line component (km s"!); (g) 
uncertainty on the velocity width (km s" l); (h) integrated line flux (10" erg cm" ^); (i) error on the line 
flux (10-15 erg s-i). 



Line 


Component 


Rest 


Velocity 


AVelocity 


Velocity 


AFWHM 


Line 


AFlux 






Wavelength 


shift 


shift 


Width 




Flux 








(A) 


(km s-i) 


(km s-i) 


(km s-i) 


(km s" 1) 


(10-15) 


(10-15) 


a) 


(6) 


(c) 


(d) 


(e) 


{/) 


(s) 


(h) 


(i) 


[O III] 


n 


5006.9 








340 


23 


1.37 


0.09 




i 


5006.9 


-402 


9 


1255 


12 


19.78 


0.27 




b 


5006.9 


-1980 


36 


1944 


65 


7.97 


0.31 


[OI] 


n 


6300.3 








340 




0.52 


0.03 




i 


6300.3 


-48 


9 


1124 


33 


2.32 


0.10 




b 


6300.3 


-1095 


66 


2671 


91 


2.37 


0.13 


[S II] 


n 


6716.4 








340 




1.28 


0.03 




i 


6716.4 


-222 


90 


1124 


119 


1.42 


0.34 




b 


6716.4 


-1374 


213 


2857 


229 


0.79 
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structure. This, along with the strong asymmetry in spa- 
tial extent of the emission, suggests the ambient gas has not 
settled down into a relaxed state. Figure 3 also reveals large 
velocity shifts (see section 4.2), up to ~ 2000 km s~^, be- 
tween the different components fitted in the nuclear aperture 
- shifts too large to be explained by gravitational motions 
and are likely to result from jct-cloud interactions, quasar 
induced winds or starburst-drivon superwinds. 

Wc assume the narrowest component in the nucleus 
traces the rest frame of the galaxy for five reasons: 

(i) The narrow component is consistent with the narrow 
component in the extended [O 11] emission, the only compo- 
nent which can be traced across the entire spatial range. 

(ii) The rcdshift of the narrow component is consistent 
with the deep HI 21cm absorption observed by Mirabel 
(1989) and Morganti et al. (2002) which is assumed to em- 
anate from the ambient, quiescent ISM in the galaxy. It is 
also consistent with the redshift of the stellar absorption 
lines (Grandi 1977). 

(iii) The narrow component is by far the least kinemati- 
cally disturbed component (FWHM = 340 ± 23 km s"^). 

(iv) The narrow component has a small velocity ampli- 
tude (Ave! < 250 km s^i) consistent with gravitational mo- 
tions in low-z galaxies (Tadhunter, Fosbury & Quinn 1989). 

(v) As will be shown in the next few sections, whilst dif- 
ferent models for the intermediate and broad components 
are required to reproduce the profiles of different emission 
lines, the narrow components remain the same for all mod- 
els and all emission lines. The narrow component data are 
presented in Table 2. 

4.2 Modelling the emission lines 

In section 4.1 wc showed the narrow component of [O II] 
traces the rest frame of the galaxy. Here we model all strong 
emission lines in the nucleus to search for outflows. 

Prior to the modelling of the emission line profiles in 



the nuclear regions a continuum model, comprising an el- 
liptical galaxy template, a power-law and a nebular contin- 
uum, was subtracted from the data. The emission lines were 
modelled using the minimum number of Gaussians required 
to produce a physically viable good fit. In order to model 
the unusually broad emission lines in the nucleus of PKS 
1345+12, it is essential to use three Gaussians (note by nu- 
cleus we mean an aperture 2.16 arcsec wide centred on the 
continuum emission). Figure 4 shows the best fit model for 
the [O 11I]AA4959,5007 lines. 

The [O 111] lines were fitted using three constraints 
in accordance with atomic physics: the flux ratio between 
[O III]A4959 and [O III]A5007 was set at 2.99:1 (based on 
the transition probabilities); the widths of the correspond- 
ing components of each line were forced to be equal; and 
the shifts between the corresponding components of each 
line were flxcd to be 48. OA. (Note, the fitting program used 
can only work with a wavelength difference and not a ra- 
tio of wavelengths. However, we find the incurred error is 
smaller than our estimated uncertainty.) The best-fltting 
model comprises 3 components for each line: 

(i) A narrow component, FWHM = 340 ± 23 km s"^ 

(ii) An intermediate component, FWHM = 1255 ± 12 
km s^^, blue shifted by 402 ± 9 km s^^ with respect to the 
narrow component. 

(iii) A broad component, FWHM = 1944 ± 65 km s"\ 
blue shifted by 1980 ± 36 km s~^ with respect to the narrow 
component. 

AH three components are plotted in Figure 3. 

Initially, in modelling the other emission lines, we as- 
sumed that one model would reproduce all emission lines, a 
technique which has been successful in studies of jet-cloud 
interactions in powerful radio galaxies (e.g. Villar-Martin 
ct al. 1999a, Villar-Martin et al. 1999b). Hence, we at- 
tempted to model the other emission lines with the same 
three components and the same velocity widths and shifts 
as [O III], leaving the relative fluxes in the kinematic sub- 
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Figure 4. Model for [O ill] AA4959,5007 in the nucleus. The best 

fitting overall profile is ovcrplottcd (bold line) on the extracted 
spectrum (feint line). All of the 6 components (3 for each line) arc 
also plotted: the dot-dashed lines trace the narrow components, 
the solid lines represent the intermediate components and the 
dotted lines trace the broad components. 



components to vary. We will refer to this as 'the [O III] 
model' hereafter. 

The [O III] model gives a good fit to H/3 and is there- 
fore assumed to reproduce Ha. However, as demonstrated 
in Figures 5 and 6, the [O III] model poorly reproduces the 
emission line profiles for [O I]AA6300,6363 and, to a lesser 
extent, [S II]AA6716,6731. 

For [O I]AA6300,6363 and [S II]AA6716,6731, various 
models were tried ranging from free fitting (i.e. only the shift 
between the two [O I] or [S II] lines, the requirement that 
the corresponding components had the same width and, in 
the case of [O I], the intensity ratio 6300/6363 = 3.00, were 
fixed) to highly constrained models to find the best overall 
model. For all models, the narrow component varied little 
and was consistent in velocity width and shift with the nar- 
row component of the [O III] model. Hence, for all models 
for all lines, the narrow component width and shift were 
constrained to be the same as that in the [O HI] model. 
Note that, it is essential to use 3 Gaussians to reproduce all 
emission line profiles, irrespective of the model used. 

Figure 5 shows the best fit model for [O I]AA6300,6363 
(the [O I] model). Apart from constraints for the narrow 
component, only the shift between the two lines (63.5 A) 
and the intensity ratio (A6300/A6363 = 3.00), known from 
atomic physics, were constrained. The [O HI] model for [O 
I] is also overplotted and clearly does not reproduce the 
emission lino profile well. 

Our model for [O I] assumes that no emission from [S 
III]A6312 and [Fc X]A6376 is significant. Models for [O I] in- 
cluding [S III]A6312 and [Fe X]A6376 were tried with the [O 
HI] model for [O I]. Whilst the overall line profile was mod- 
elled well, the contribution of [S III]A6312 and [Fe X]A6376 
became unreasonably large and the model was discounted. 

Similarly, Figure 6 shows the best fit model for the den- 
sity diagnostic [S II]AA6716,6731. Again, the [O HI] model 



does not provide a good fit, but the discrepancy is not as 
extreme as for [O I] AA6300,6363. 

Due to the complex, broad emission lines in the nucleus 
and uncertainties in the continuum subtraction, it was diffi- 
cult to fit the fainter lines in the nucleus and so the results 
are omitted from this paper. 



4.3 Electron densities 



During the early stages of radio source evolution, the 
nuclear regions of radio sources are likely to harbour a rich 
ISM, possibly left over from a recent merger. Here we in- 
vestigate the density of the line emitting gas using [S II] 
6716/6731, the density diagnostic ratio. 

In section 4.2 and Figure 6 we showed that the [O HI] 
model fit to the [S II]AA6716,6731 doublet is not wholly sat- 
isfactory. The initial free fit for [S II] (i.e. only the narrow 
component width and position was constrained) produced 
a good fit to the overall line profile, but the ratio was out- 
side the range of allowed values for both the intermediate 
and broad components ([S II] 6716/6731 < 0.44 - the ratio 
at the high density limit). For the next attempt, the same 
constraints were applied with the addition of the ratios of 
the broad and intermediate components set at the high den- 
sity limit (6716/6731 = 0.44). This produced a good overall 
fit for the doublet which was physically viable (the narrow 
component ratio was well within the allowed range). 

For subsequent fits, the [S II] 6716/6731 intensity rar 
tio for the broad and intermediate components was slowly 
increased (i.e. the density decreased) until the model broke 
down and no longer fitted the blend well. Hence we esti- 
mate upper limits for the intensity ratio [S II] 6716/6731 of 
0.57 and 0.60 for the intermediate and broad components re- 
spectively. Reassuringly, the density estimates obtained by 
fitting the [S II] doublet with the [O HI] model were similar 
to those obtained using the [S II] model. 

The density was calculated using the IRAF program 
TEMDEN which is based on the five-level atom calculator 
developed by dc Robertis, Dufour & Hunt (1987), assuming 
an electron temperature of 10 000 K. Hence, upper limits 
for the intensity ratio [S II] 6716/6731 correspond to lower 
limits in density and we calculate densities of rie > 5300 
cm~^ and rie > 4200 cm~^ for the intermediate and broad 
components respectively. The [S II] 6716/6731 ratio for the 
narrow component is consistent with the low density limit 
and we estimate an upper limit of Ue < 150 cm~^ for the 
narrow component. 



4.4 Reddening 

In section 4.3 we showed that the nuclear regions of PKS 
1345+12 harbour a dense ISM (ue > 4000 cm"^). Observa- 
tionally, large amounts of gas and dust will be detected as 
reddening of the optical spectrum. 

We have investigated reddening in PKS 1345+12 us- 
ing three independent techniques: the Ha/H/J ratio, the 
Paa/H/3 ratio and nebular continuum subtraction, assuming 
a simple foreground screen model for interstellar dust. 
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Figure 5. Models for [O I] AA6300,6363 in the nucleus. The bold line represents the best fitting 
model for [O I] (the [O I] model) and the three components for each line are plotted (see Figure 4 
for key). The dashed line represents the [O III] model for [O I] (see section 4.2). 
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Figure 6. Models for [S II]AA6716,6731 in the nucleus. The bold line represents the best fitting 

model for [S II] (the [S II] model) and the three components for each line arc plotted (see Figure 4 
for key). The dashed line represents the [O III] model for [S II] (see section 4.2). 
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Figure 7. The nebular continuum, a) The bold line traces a section of the original nuclear spectrum of PKS 1345+12 in the blue. The 
dotted line represents the total nebular continuum (pseudo + theoretical case B continua for all three components) created using the 
measured H/3 flux (i.e. assuming zero reddening). The feint line traces the nebular-continuum subtracted spectrum. There is clearly a 
sharp discontinuity at ~ 3650 A in the nebular continuum subtracted spectrum corresponding to the jump in the nebular continuum 
at the same point. This indicates that the measured values of HyS are highly reddened resulting in an incorrectly calculated nebular 
continuum, b) The same key as for a) but the nebular continuum subtracted was generated using the reddening corrected B./3 fluxes 
from the E(B-V) values calculated from the Ha/H/3 ratios and then rc-rcddcncd accordingly (sec section 4.4 and Table 3). This confirms 
the E(B-V) values for reddening are of the correct order as there is no longer a sharp discontinuity at ~ 3650 A. 



4-4- i Balmer line ratios 

The degree of reddening can be estimated using the Balmer 
line ratio Hq/H/?. Note that the Ha/[N II]AA6548,6583 
blend in the nuclear spectrum of PKS 1345+12 is highly 
complex - all three lines are broad, have blue asymmetries 
and require three components each. Initially, both the Hq 
and [N II] lines were fitted using the [O III] model. The [O 
III] model was used for Ha since H/3 is well modelled by it. 
For [N II], various models were tried and the [O III] model 
produced the best overall fit for the blend. Using this [O III] 
model we measured Ha/H/9 ratios of 3.32 ± 0.33, 5.26 ± 
0.28 and 18.81 ± 4.74 leading to E(B-V) values of 0.06 ± 
0.05, 0.42 ± 0.10 and 1.44 ± 0.50 for the narrow, intermedi- 
ate and broad components respectively (see Table 3; Seaton 
1979 extinction law assumed). 

Given the complexity of the blend, we also re-modelled 
the lines to find the lowest possible flux for Ha, hence a 



lower limit for the ratio Ha/H/3 and the E(B-V) value. We 

used the [O HI] model for Ha but for ]N II] only constrained 
the narrow components and for the intermediate and broad 
components, the flux ratio ([N II] 6583/6548 = 3.00) and 
the shift between the lines (35.27 A) and forced the velocity 
widths of the corresponding components to be equal. This 
technique gave lower limits for Ha/H/3 of Ha/H/3 > 3.91 
and Ha/H/3 > 9.81 for the intermediate and broad compo- 
nents respectively corresponding to lower limits for E(B-V) 
of E(B-V) > 0.19 and E(B-V) > 0.92 for the intermediate 
and broad components respectively. 



4-4-^ Comparison of the optical and infra-red 

Vcilloux, Sanders & Kim (1997) detected two components to 
Paa: broad (FWHM - 2600 km s"^) and narrow (FWHM 
~ 900 km s^'^), and they identified the broad component as 
broad line region (BLR) emission from a lightly-extinguished 
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Table 3. Reddening values. Columns are: (a) Quantity, where: Ha/H/3 is the flux ra- 
tio between the measured values of Ha and H|S; E(B-V)Ha/H/3 is the E(B-V) value 
calculated from Ho/H/3 and the standard interstellar extinction curve (Seaton 1979); 
E(B-V)PaQ/H,3 is the E(B-V) value calculated from the ratio between the value of Paa 

quoted by Veillcux, Sanders & Kim (1997) and the measured H/3 flux; Ffj^ is the red- 
dening corrected (assuming the E(B-V) value calculated using the Ha/H/3 ratio) H/3 line 
flux (10^"'"''' erg s^-' cm^^); Ljj/j is the rest frame H/3 luminosity (10^" erg s^^) and Mgas 
is the mass of line emitting gas in solar masses; (6) value of particular quantity for the 
narrow component; (c) error for narrow component; column pairs (d) and (e) and (/) 
and (g) are as for (6) and (c). *: using quoted Paa flux from Veilleux, Sanders & Kim 
(1997). 
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quasar nucleus. Although this interpretation is consistent 
with the detection of a point source in high-resolution near- 

IR observations (Evans ct al. 1999), there must be some 
doubt about it given the detection of broad components to 
the forbidden lines at optical wavelengths. Unfortunately, it 
is difficult to establish the link between the various oirtical 
and infrared kinematic components with greater certainty, 
because Veilleux, Sanders & Kim (1997) neither publish the 
line shifts, nor fit an intermediate width component. How- 
ever, by assuming that all the flux in the broad Paa is as- 
sociated with the same component that emits the broad H/3 
detected in our optical spectrum and by varying this flux 
by a factor of 2 to account for possible slit-loss differences 
between the observations, we obtain an upper limit for red- 
dening of the optical broad component of E(B-V) < 2.0 - 
consistent with the other estimates obtained above. 



4.4-3 The nebular continuum 

Another indication of large reddening in the nuclear spec- 
trum of PKS 1345-1-12 comes from consideration of the neb- 
ular continuum, calculated as part of the continuum mod- 
elling process prior to measurement of the emission lines. 
The nebular continuum for each kinematic component was 
created using the technique described in Dickson et al. 
(1995). The nebular contirmum strength is not particularly 
sensitive to temperature but we have assumed an electron 
temperature, Te = 10 000 K. 

Initially, we assumed zero reddening for all components. 
Figure 7a) shows the original spectrum, the zero reddening 
nebular continuum and the nebular continuum subtracted 
spectrum in the nucleus of PKS 1345-1-12 in the blue. After 
subtracting the nebular continuum, a clear discontinuity is 
left at the position of the Balmer edge at ~ 3645A which 
is unphysical. This indicates that the nebular contirmum 
has been over-estimated in the zero reddening case and that 
reddening is likely to be significant in the nucleus of PKS 
1345-1-12. 



Using the E(B-V) values calculated from the measured 
Ha/H/3 ratios (see section 4.4.1) and the standard inter- 
stellar extinction curve (Seaton 1979), we corrected the H/3 
fluxes for reddening and calculated the reddened nebular 
continuum spectra. 

Figure 7b) shows the original spectrum, the new neb- 
ular continuum (assuming E(B-V) values of 0.06, 0.42 and 
1.44 for the narrow, intermediate and broad components re- 
spectively) and the nebular continuum subtracted spectrum 
in the nucleus of PKS 1345-1-12. Clearly, when the reddened 
nebular continuum is subtracted, there is no longer a discon- 
tinuity. This is further evidence for large reddening in the 
nucleus of PKS 1345+12 

To further investigate the degree of reddening in the 
broad and intermediate components (the amount of redden- 
ing in the narrow component is negligible) , this process was 
repeated for various E(B-V) values for the intermediate and 
broad components. The intermediate component contributes 
most of the flux in the nebular continuum and we estimate 
a lower limit of E(B-V) > 0.3 for this component with some 
confidence - at this level of reddening, the discontinuity in 
the nebular continuum subtracted spectrum disappears. 

The variation of the broad component had negligible 
effect on the resultant nebular continuum. Hence, the lower 
limit of E(B-V) > 0.92 for the broad component comes from 
the calculated Balmer line ratios. 

To summarise, the nuclear regions of PKS 1345-1-12 har- 
bour a rich ISM and are subject to large reddening. Each of 
the three independent techniques for estimating reddening 
give consistent results and it is clear there is a strong corre- 
lation between the broadness of component and the degree 
of reddening, in the sense that the broadest components are 
the most highly reddened. 



4.5 Estimating the gas mass 

We have seen that the nuclear regions of PKS 1345-1-12 har- 
bour a rich, dense ISM, but is there sufficient mass to confine 
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and frustrate the radio source? Emission line luminosities 
are related to the mass of line emitting gas by: 

L(H/3) 



M 



where Ne is the electron density (cm~^); mp is the mass of 
the proton (kg); L(H/?) is the luminosity of the H/3 line (erg 
is the effective recombination coefficient for H/3 



e// 



(cm^ s~^) and b.UH0 is the energy of an H/3 photon (erg). 

In section 4.3 we calculated densities of rie < 150 cm~^, 
ne > 5300 cm^"^ and Ue > 4200 cm~^ in the regions emit- 
ting the narrow, intermediate and broad components respec- 
tively. Using the presented equation, lower limits in density 
correspond to upper limits in mass. Hence, we estimate lim- 
its on the mass of line emitting gas of Mgas > 2.61 x 10* 
Mo, Mgas < 0.92 X lO'' Mo and Mgas < 5.64 x lO'"^ Mo 
in the regions emitting the narrow, intermediate and broad 
components respectively (see Table 3). This leads to an up- 
per limit of order 10® M0 for the total mass of line emitting 
gas in the kinematically disturbed emission line components. 



5 DISCUSSION 

5.1 A model for PKS 1345-M2 

The broad, multiple-component emission lines in PKS 
1345-1-12 indicate that line emission emanates from at least 
three kinematically distinct regions in the NLR along the 
line of sight. To determine what these components repre- 
sent, it is necessary to establish the exact redshift of the 
galaxy rest frame. 

In section 4.1, we argued that the narrow component 
represents the ambient, quiescent ISM in the galaxy halo. 
We are confident of this result - the narrow component is 
the only component consistently found in all emission lines, 
irrespective of model required or spatial position; kinemati- 
cally, it is the least disturbed component and it is consistent 
with the deep, narrow HI 21cm absorption trough (Mirabel 
1989, Morganti et al. 2002) usually associated with the am- 
bient, quiescent gas halo surrounding radio sources. Hence, 
the blue asymmetries in the emission line profiles represent 
material moving towards us, blue shifted with respect to the 
galaxy. 

In sections 4.3 and 4.4 we showed that the intermediate 
and broad components originate from regions that are dense 
(ue > 4000 cm"^) and dusty (E(B-V) > 0.30 and E(B-V) 
> 0.92 for the intermediate and broad components respec- 
tively). It is unlikely, therefore, that we would see emission 
from material on the far side of the nucleus - the further 
away from the observer in PKS 1345-1-12 along the line of 
sight, the more extinction the emission will suffer. The blue 
shifted material in PKS 1345-1-12 is therefore likely to be 
material in outflow on the side of the nucleus closest to the 
observer. 

Indeed, similar conclusions have been made about the 
compact, flat-spectrum radio source PKS 1549-79 by Tad- 
hunter et al. (2001). The optical spectrum of PKS 1549-79 
shows many characteristics similar to PKS 1345+12, namely 
blueshifted broad emission lines. Moreover, a large extinc- 
tion for this source is suggested by the fact that, whilst the 
flat radio spectrum indicates the presence of a QSO, there is 



no evidence at optical wavelengths for the broad permitted 
lines and non-stellar continuum characteristic of QSOs. We 
will now build on the simple model proposed for PKS 1549- 
79 to account for more complex characteristics observed in 
PKS 1345+12. 

In sections 4.2 and 4.3 we showed that three Gaussians 
were required to fit all the emission lines, irrespective of the 
model required to reproduce the profile. The models vary in 
both velocity width and shift of the intermediate and broad 
components with respect to the narrow component. We in- 
terpret this as three kinematically and physically distinct 
regions in PKS 1345+12, each responsible for the emission 
of one of the components (i.e. narrow, intermediate or broad) 
although gradients in ionisation/density, coupled with gra- 
dients in velocity must exist across these regions to produce 
different broad or intermediate components (see Figure 8). 

We can determine the locations of the three kinemati- 
cally distinct regions using the reddening results. In section 
4.4 we showed that whilst there is significant reddening in 
the nuclear spectrum of PKS 1345+12, the amount of red- 
dening in the three components differs substantially. The 
narrow component has negligble reddening and must there- 
fore originate from the outer regions of the galaxy. This is 
consistent with other evidence suggesting the narrow com- 
ponent originates from the galaxy halo (see section 4.1 and 
above). Similarly, the intermediate and broad components 
are subject to significantly more reddening with the degree 
of reddening correlating with the broadness of the compo- 
nent. The broadest component is therefore likely to origi- 
nate from the regions closest to the galactic centre and the 
intermediate component slightly further out but within the 
narrow component emitting region. 

We now propose an evolutionary scenario for PKS 
1345+12: 

(i) PKS 1345+12 underwent a merger event in its recent 
history. This merger consisted of two galaxies (PKS 1345+12 
has a double nucleus), at legist one of which was gas rich. The 
merger event started so recently it is still ongoing, indicated 
by the presence of two nuclei (e.g. Gilmore & Shaw 1986; 
Heckman et al. 1986; O'Dea et al. 2000). This is consistent 
with the idea that many radio galaxies appear to have been 
triggered by merger events. 

(ii) The merger event provided large amounts of gas and 
dust to the regions surrounding the nucleus, obscuring it 
from view. 

(iii) The current activity was triggered in the the western 
nucleus, the nucleus coincident with the radio source. Note, 
we use the term 'current activity' because, whilst we can say 
with confidence that the radio source is young (see section 
5.2), we cannot rule out the possibility of previous activity 
phases which may be responsible for the extended, diffuse 
halo emission. 

(iv) While the active nucleus is obscured, the outer parts 
(halo etc) of the galaxy start to settle down. The halo is re- 
sponsible for the emission of the narrow component and the 
deepest HI absorption trough (Mirabel 1989; Morganti et al. 
2002). Hence, the velocity width of the narrowest component 
is consistent with gravitational motions. PKS 1345+12 is a 
GPS source and is therefore observed relatively shortly af- 
ter the triggering of the radio activity. With time, the radio 
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Figure 8. Schematic diagram of the model for PKS 1345+12. 



source will expand, first becoming a larger compact steep 
spectrum radio source (CSS) before becoming extended. 

(v) As the radio jets expand, material is shocked, accel- 
erated and pushed aside as the jets hollow out cavities on 
either side of the quasar (bi-polar). As well as jet-cloud inter- 
actions, quasar- or starburst-induced winds arc also likely to 
be important at this stage. Indeed, outflow components with 
velocity shifts similar to those measured in PKS 1345-1-12 
have also been detected in ULIRGs without strong radio 
emission (Wilman, Crawford & Abraham 1999). Therefore, 
the dominant driving mechanism for the outflows is cur- 
rently unknown. 

(vi) Eventually, the radio jets will tunnel through the co- 
coon and expand freely to evolve into a large scale, extended 
radio galaxy. 



5.2 PKS 1345-1-12 - young radio source 

A crucial question regarding compact radio sources involves 
the amount of material in their circumnuclcar regions. In 
the evolutionary scenario outlined above, we have assumed 
that PKS 1345-1-12 is a young radio galaxy. It is clear that 
the circumnuclear ISM in PKS 1345-1-12 is both dense and 



dusty, but is there sufficient mass to confine and frustrate 
the radio source? 

In section 4.5, we estimated the mass of line emitting 
gas using the reddening corrected H/3 luminosities (see Tar 
ble 3). We estimated an upper limit for the total mass of 
line emitting gas (i.e. for the kinematically disturbed com- 
ponents combined) of < 10^ M©. 

Several simulations have been run varying the distribu- 
tion and density of the instorstollar medium to determine the 
mass of gas required to frustrate a radio source. Hydrody- 
namical simulations carried out by De Young (1993) showed 
that for low to intermediate power radio sources (10*^ - 10'''' 
erg s~^) in both a homogeneous and a two-phase ISM, gas 
masses of order lO" Mq wore sufficient to confine the ra- 
dio source for its entire lifetime. For the high power radio 
sources (10*^ erg s~^), the masses required would be too 
large to be credible. Analytical models of jet interactions by 
Carvalho (1994,1998) suggest confinement is possible with 
much lower eas masses, of order lO" - 10^" M., if the ISM 
is clumpy and within a scale of hundreds of pc of the radio 
source. 

At face value, these sirrmlations, considering cither a 
smooth medium or small clouds evenly distributed through- 
out a sphere, may lead us to conclude that our estimate of 
Mgas < 10® Mq may be too small to confine and frustrate 
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the radio source and that PKS 1345+12 is a young radio 
source. However, if the ISM forms large dense clouds, much 
smaller masses could significantly disrupt the propagation of 
the jet (e.g. forming bubbles as in M87; Bicknell & Begelman 
1995) which may lead to some degree of frustration. 



6 CONCLUSIONS AND FUTURE WORK 

PKS 1345+12 is most likely a young radio source observed 
sometime after the start of the merger event which triggered 
the nuclear activity. We observe complex, multiple compo- 
nent emission lines in the nucleus which, due to largo red- 
dening in the nucleus, we interpret as material in outflow. 
Line emission originates from three kinematically distinct 
regions over which exist gradients in velocity, acceleration 
and possibly ionisation potential. The major outstanding 
issues for PKS 1345+12 concern the acceleration and ion- 
isation mechanisms for the gas. In particular, whether the 
outflows are driven by the jets, AGN driven winds or star- 
burst driven super-winds. In future, it should be possible 
to address these issues using high resolution imaging obser- 
vations of the outflow regions in combination with detailed 
studies of the emission line ratios derived from our spectra. 

Overall, our results provide evidence that the activity 
has a major impact on the circum-nuclear gas in the early 
stages of radio source evolution. 
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